D-wave bosonic pair in an optical lattice 



O 
(N 



O 
(N 



C/3 



I 

o 



(N 
> 
00 

O 
in 
o 



o 



X 



Zi Cai\ Lei Wang^ , Jian Li^ Shu Chen^ X. C. Xie^'^.s^ ^ Yupeng Wang^ 
^Department of Physics, University of California, San Diego, California 92093, USA 
^ Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100080, P. R. China 
Texas Center for Superconductivity and Department of Physics, 
University of Houston, Houston, Texas 77204, USA 
'^International Center for Quantum Materials, Peking University, Beijing 100871, China and 
^Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078, USA 

(Dated: March 22, 2011) 

We present a bosonic model, in which two bosons may form a bound pair with d-wave symmetry 
via the four-site ring exchange interaction. A d-wave pairing superfiuid as well as a d-wave density 
wave (DDW) state, are proposed to be achievable in this system. By the mean field approach, we 
find that at low densities, the d-wave pairs may condensate, leading to a d-wave bosonic paired 
superfiuid. At half filling, a d-wave Mott insulator could be realized in a superlattice structure. At 
some particular filling factors, there exists a novel phase: d-wave density wave state, which preserves 
the d-wave symmetry within plaquette while spontaneously breaks the translational symmetry. The 
DDW state and its corresponding quantum phase transition in a two-leg ladder are studied by the 
time-evolving block decimation (TEBD) method. We show that this exotic bosonic system can 
be realized in the EEC zone of cold Fermi gases loaded in a two-dimensional (2D) spin-dependent 
optical lattice. 

PACS numbers: 05.30.Jp, 03.75.Nt, 74.20.Mn, 73.43.Nq 



Recently, ultracold atoms in optical lattice have pro- 
vided a perfect platform for simulating quantum many- 
body model in condensed matter physics. Because of 
the flexible tunability of parameters such as the hop- 
ping amplitudes, interaction or even the dimensionality 
of the system, ultracold atomic systems allow us to di- 
rectly study some fundamental Hamiltonian systems and 
their associated phase transitions, such as boson Hub- 
bard model and the superfluid to Mott insulator phase 
transitioi*^, or the recent realization of the repulsive 
or attractive fermionic Hubbard model^. In addition, 
some exotic phases emerged from the low-dimensional 
strongly correlated systems, such as resonating valence 
bond (RVB) stated"—, d-wave superfluidity^-'', deconfined 
Coulomb phase^iS as well as the topological insulator 
with fractional statistic and topological orderi^^iScan also 
be investigated in the cold atom systems. Further- 
more, the uniqueness of cold atomic system also pro- 
vides new playgrounds for physicists, such as the strongly 
correlated model for higher spin systems, higher orbital 
systemsiir— or for the optical superlattice^ii^"— . A two- 
dimensional (2D) optical superlattice may be constructed 
by imposing two optical lattices with different periods to 
form an array of plaquettes. The hopping amplitude and 
interaction for atoms between these plaquettes are much 
smaller than that within the plaquette. One of the ex- 
otic phase emerges in the superlattice is the d-wave Mott 
insulatori^iii and d-wave superfluid^®. The d-wave Mott 
insulator is the insulator state with local d-wave symme- 
try, i.e., if we rotate the site within a plaquette by tt/2, 
the wave function reverses its sign. When we introduce 
holes into the d-wave Mott insulator, two holes tend to 
bind together within the plaquette to form a Cooper pair 
with local d-wave symmetry, and the propagation of the 



d-wave pairs between different plaques leads to the d- 
wave superfluid. 

The mechanism of pairing with d-wave symmetry has 
played an important role in the high-Tc superconduc- 
tor. Though without rigorous proof, numerous evidences 
strongly support the existence of d-wave superconductor 
(or superfluidity in cold atomic system) near the half- 
filling. The background of Neel state with antiferromag- 
netic correlations plays a key role in this mechanism of 
d-wave symmetrjsi^. In cold fermionic atom system, how- 
ever, the binding energy of a d-wave pair is much smaller 
than the hopping amplitude, which makes it difficult to 
directly simulate the d-wave mechanism of the high-Tc 
superconductor in cold atomic systems. One solution is 
to trap the d-wave pair within a plaquette of the optical 
superlattice. In this paper, we propose a novel mecha- 
nism for realizing the d-wave pairing. Different from that 
in high-Tc superconductor induced by antiferromagnetic 
correlation, the d-wave paring here is induced by the four- 
site ring exchange interaction. 

Before discussing the physical realization of the effec- 
tive Hamiltonian, we first present the Hamiltonian, which 
is a hard-core bosonic model with a strong nearest neigh- 
bor (NN) repulsive interaction and a four-site ring inter- 
action: 



H = [tal aj + VniTij] — fi 
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where {ij) denotes a pair of nearest-neighbor sites and 
{ijkl) are sites on the corners of a plaquette. In this pa- 
per ,we focus on the parameters region V K ^ \t\ > 0. 
The four-site ring exchange interaction with positive co- 
efficient (K>0) may be very important in determining 
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FIG. 1: (a) Anisotropic hopping matrix elements in the 2D 
anisotropic spin-dependent optical lattice (green t; red l): tjv 
and ta- (b) The effective hopping term for the hardcore boson 
as a result of second order perturbation, (c) The effective 
interaction term for the hardcore boson, (d) The four-site 
ring exchange interaction for hardcore boson as a result of 
fourth perturbation, (e) Two eigenstate s-wave and d-wave 
symmetry within one plaquette 



the properties of this system (at least at low densities) 
and leads to the exotic boson pairs with the d-wave sym- 
metry. Though in most systems, the effect of this four- 
site ring exchange interaction is much smaller compar- 
ing to the two-site hopping amplitudes because it usu- 
ally comes from the fourth order perturbation, below we 
would show that we can realize the Hamiltonian as well as 
the corresponding parameter region (y ^ if ^ |i| > 0) 
in the BEC zone of cold Fermi gases loaded in a two- 
dimensional (2D) spin-dependent optical. Similar model 
without the NN repulsive interaction has been proposed 
to study the exotic phases in cold atom system such as the 
deconfined phased or Bose-metal phased"—. However, as 
we will show below, the strong NN repulsive interaction 
in Eq.([T]) makes the novel Bose-metal phase unstable. 

First, we would discuss the experimental realization 
of our Hamiltonian. (1) as well as the corresponding pa- 
rameter regions. Most of above discussion are based on 
the parameter region: y ^ A' ^ |i| > in Hamilto- 
nian. ([T]). However, in most systems, no matter in solid 
physics or cold atom physics, the effect of this four-site 
ring exchange interaction is much smaller comparing to 
the two-site hopping amplitudes because it usually comes 
from the fourth order perturbation. Below we would 
show that not only the Hamiltonian. ([1]) but also the cor- 
responding parameter region (y ^ X ^ |f| > 0) could 
be realized in the BEC zone of cold Fermi gases loaded 
in a two-dimensional (2D) spin-dependent optical^. A 
fermionic Hubbard model with spin- dependent hopping 
has been proposed^i, by tuning the lasers between hype- 
fine structure levels of ^'^-R' atoms. The tunneling matrix 
elements for the two spin components are spin-dependent 
and can be tuned with different anisotropy. In our case, 
we tune the hopping so that spin up | f) atoms prefer 



to hop along the x axis and spin down | J,) atoms pre- 
fer to hop along the y axis, which means there are two 
kind of typical hopping amplitude ta and i^. ta repre- 
sents the hopping amplitude for the | f) (| |)) fermions 
along X (y) axis, while t-^ denotes the hopping amplitude 
for the I t) (I i)) fermions along y (x) axis (as shown 
in FigUJa)). The ratio 5 = tT^jta can be tuned experi- 
mentally and we choose the high anisotropic condition: 

<ti ta- In addition, we can use Feshbach resonances to 
manipulate the interactions and adjust it from repulsive 
to attractive. The unconventional pairing with attrac- 
tive interaction in the anisotropic spin-dependent optical 
lattice has been analyzed recently^^. 

We load the fermions into the 2D spin-dependent op- 
tical lattice defined above. Then we use Feshbach reso- 
nances to make the two fermions occupying the same sites 
binding together to form a bosonic molecule. Obviously 
this molecule is hard-core in nature. We assume that 
the binding energy is large enough that all fermions are 
tightly bound into bosonic molecule and the system en- 
ter a BEC zone. Next we will analyze the dynamics and 
interactions of these new bosons to show how can we con- 
struct the Hamiltonian ([!]) as well as the corresponding 
parameter region in this spin-dependent optical lattice. 

As shown in Fig|TJ both the two-site hopping and in- 
teraction involve the second order perturbation via a vir- 
tual process. Taking the hopping term for example (Fig.l 
(b)), from the standard second order perturbation theory 
we can obtain the effective hopping amplitude of the bo- 
son: t = —tatTr/U, ta and tTr has been defined above. U is 
the binding energy for a bosonic molecule formed by two 
fermions. Similarly, the virtual process in Fig.l (c) plays 
a role similar to the term VSfSj {V = t^/U > 0) in the 
spin model, which means the NN repulsive interaction 
VniJij in the boson language (those terms proportional 
to rii are absorbed into the chemical potential). Because 
of the strong anisotropic hopping of the fermions in the 
spin-dependent optical lattice: t^r "^ta, we have V ^ t. 

The virtual process shown in Fig. 1(d) is the leading 
term from the fourth order perturbation. Apparently it 
involves four sites within a plaquette and thus results in 
a ring exchange interaction in Eq[TJ From the standard 
perturbation theory, we can get K — t*/U^ > 0, which is 
much larger than the contribution of all the other four- 
site virtual processes. Due to the strongly anisotropic 
hopping of the fermions in the spin-dependent lattice, it 
is possible to adjust the parameters (6 — t-^/ta — > 0) in 
Hamiltonian ([T]) to satisfy K ^ t , which means that 
the pair binding energy is much larger than the single- 
particle hopping energy. 

To define the symmetry of a pair, we introduce a lo- 
cal operator Dp, which rotates the four sites within the 
plaquette p cyclically by an angle 7r/2. The sign of K is 
important for the symmetry of this bosonic pair. When 
if < 0, it is s-wave {Dp — 1); while if > is d-wave 
(Dp = —1). This can be seen with just one plaquette. 
Loading two bosons into one plaquette: due to the strong 
NN repulsive interaction, the only two possible configura- 
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FIG. 2: The mean-field result at zero temperature (we set 
K=l and t=0.1). The dependence of (a) energy gap A and 
(b) the energy of the ground state Eg of the d-wave pairing 
phase (solid line), and single particle BEG (dashed line) on 
the total particle number. 



tions in the plaquette is two bosons occupying the diago- 
nal sites. There are two eigenstates of the ring-exchange 
term, denoted as \s) and |(i)(Fig[Ue)), with eigenvalues 
K and -K, respectively. We neglect the single particle 
hopping term because K ^ t. Notice that when K > Q, 
the ground state of this plaquette is a d-wave state . At 
low densities, since the pair binding energy is much larger 
than the single particle hopping energy, the bosons prefer 
to move as pairs rather than hopping independently. At 
low density, these bosonic pairs with d-wave symmetry 
will condensate to form a d-wave superfluid. 

A natural question arose here is whether the ring ex- 
change interaction would make the system to form bound 
state with more than two bosons? Without the NN repul- 
sive interaction, this is true: at low density the ring ex- 
change, just like an attractive potential, causes not only 
two but many bosons to clump together, which leads to 
phase separations^ between an isolated boson metal clus- 
ters and vacuum. In our case, however, the strong repul- 
sive interaction will make the isolated clusters consisted 
of more than two bosons unstable and break up into many 
bosonic pairs. 

At a low density, the interaction between the pairs is 
not important, thus, we can analyze the problem us- 
ing the mean field theory?^. We introduce a d-wave 
bosonic pair order parameter to decouple the four-site 
ring exchange interaction in our original Hamiltonian 
Eq. ([T]): {a[a\) = —{a\a}^ = A, where the minus 
sign is due to the d-wave character. As we concen- 
trate on the low-density limit, the hard-core constraint 
is expected to be irrelevant. We decouple the NN in- 
teraction by the Hartree-Fock approximation: VniUj — 
V{ni)nj -\- ni{nj) — V{ni){nj). The Hamiltonian in EqU] 
can be rewritten as: 



k 



Vn^ + h.c + 2KA'^ (2) 



Afc = 2KA sin sin ky , 
ffc = 2t(cos kx + cos ky) — fJ, + 2Vn. 



where n is the average value of the total particle number. 
The mean field Hamiltonian ^ is diagonalized by using 
the Bogoliubov transformation for bosons and we obtain 
the energy spectrum: Ek = \J {^k/'^Y ~ l^fcP- We focus 
on the zero temperature case and the ground state energy 
is given as: Eg = Y,\^Ek + 2A^K ^ [ij^-Vn-Vn^ . The 
sum is over all the k in the first Brillione zone. The 
self-consistent equations are: 



4L ^ EJk ' 



^ K -^—^ sin fcj, sin ky 

Ei, ■ 

k 



(3) 
(4) 
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FIG. 3: The optical superlattice for (a) square lattice and (b) 
ladder with periodic boundary condition 

Since we are dealing with a bosonic system, it is possi- 
ble that another BEC state with single particle conden- 
sation will compete with our d-wave pairing state. In this 
case, Eq.(3) should be replaced hy n = ric — where 
ric is the density of bosons with single particle conden- 
sation. We find that at least in our parameter regime 
t K V, there is no positive self-consistent solution 
for fic- The absence of single particle condensation has 
been observed previously in a bosonic system with cor- 
related hoppingS^. It is shown that the bosons prefer to 
pairing with each other due to the strong effective attrac- 
tive interaction. To clarify this point, we also calculate 
the energy of the single particle BEC state without d- 
wave pairing via the standard Bogoliubov approximation: 
o-i — ^/n^+Sui to decouple the Hamiltonian. (1). The re- 
sult is shown in Fig. 2 (we set K=l and t=0.1). Notice 
that at least at low densities with K ^ t, the ground 
state energy of d-wave pairing state is always lower than 
that of the single particle BEC state. 

Nextwe would turn to another limit, when the filling 
factor is 1/2. In this case, the strong NN repulsive inter- 
action induces a conventional (tt, tt) density wave phase, 
rather than the boson metal or d-wave bosonic pairs. 
However, if we load the half filling bosons into a 2D 
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FIG. 4: The dependence of (Dp) on A for a (a) 4 x 4 square su- 
perlattice (b) 2x8 ladder superlattice with periodic boundary 
condition 

optical superlattice, it is possible to recover the d-wave 
symmetry within plaquette and lead to a d-wave Mott 
insulatoriMiil. Next we will clarify this point by the exact 
diagonalization (ED) of the small size systems. We clas- 
sify all the plaquettes in the superlattice as two classes: 
P (grey plaquette in Figl3]) and P' (white plaquette) and 
the Hamiltonian in this case is given by: 

i/, = i/o + A ^ i/o, (5) 

where Hq is the Hamiltonian defined by Eq.([T]) in one 
plaquette and < A < 1. Notice that when A ^ 1, 
the situation is similar to that in a single plaquette, the 
system forms a d-wave Mott insulator with (Z?p) « — 1. 
When A « 1, the ground state should be a (tt, tt) DW with 
(Dp) = due to the strong repulsive NN interaction. 
(Dp) is the expectation value of the rotating operator 
defined above to measure the d-wave symmetry within 
one plaquette. We calculate (Dp) in a 4 x 4 superlattice 
(FigHI^a)) and a 2 x 8 ladder (FiglH^b)) with periodic 
boundary conditions, to show how it changes when we 
increase A from to 1. 

At some particular filling factor (f=l/3 for the two-leg 
ladder system and 1/4 for the 2D system), a novel phase 
emerges. The ring exchange interaction makes two bo- 
son prefer to form a d-wave pair, while the strong NN 
interaction prevents two pairs from being too close. The 
conspiracy of them makes these d-wave pairs localized 
and localized and separated as far as possible to avoid 
the strong NN interaction. The crystallization of these 
d-wave pairs leads to a novel phase: d-wave density wave 
state, which preserves the d-wave symmetry within pla- 
quette and spontaneously breaks the translational sym- 
metry, as shown in FiglSl^b) (two-leg ladder) and FiglSJc) 
(2D). Notice that unlike the half filling case, the transla- 
tional symmetry breaking is spontaneously in DDW state 
thus we don't need any superlattice structure. 

Before discussing the TEBD result, we first briefly dis- 
cuss the global phase diagram of our ladder system in 
the limit of t — in Hamiltonian. ([1]). First we analyze 
the half filling case, in the limit y — 0, it is known that 
this K-only model can be mapped to a hard-core boson 
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FIG. 5: (a) Sketch phase diagram of the model in Hamiltonian 
(1) in a two-leg ladder {t 0) (b)The structure of the DDW 
state in a two- leg ladder (1/3 filling) . (c) The structure of 
the DDW state in a 2D system (quarter-filling) 



model and could be solved exactly^'^. Its ground state is a 
gapless highly correlated state of boson: d-wave Bosc liq- 
uid (DBL). When y /C, a gapped (7r,7r) density wave 
(DW) state would dominate. A quantum phase transi- 
tion would occur when we increase the NN interaction 
V. When the fiUing factor f=l/3, in the fimit F 0, 
we anticipate that the phase of ground state is separated 
into an empty region and a half-filled region of the DBL 
state. When we increase V, as analyzed above, the com- 
petition between the NN interaction and ring exchange 
interaction would lead to a DDW state, which preserves 
the local d-wave symmetry and breaks the translational 
symmetry spontaneously. The sketch global phase dia- 
gram is shown in Fig|SJa). 

Below we would study the DDW state and the proper- 
ties of the quantum phase transition in the two-leg ladder 
system by the TEBD method. We focus on the case the 
filling factor / w 1/3. The open boundary condition is 
used to artificially shift the ground state degeneracy of 
DDW state due to the spontaneous translational sym- 
metry breaking, therefore the filling factor is not exactly 
1/3 in the ladder with finite length. For example in a 
2 X ladder, the number of boson TVf, = {2Lx. + 2)/3. 
In the thermodynamic limit, Lx ^ oo^ the filling factor 
is exactly 1/3. Because the DBL state and the phase 
separation as well as the quantum phase transition be- 
tween them have been explicitly discussed in Ref.^^, thus 
we would not discuss them here and mainly focus on 
the DDW state and the properties of the corresponding 
quantum phase transition. 

Because of the different structure of the phase sepa- 
ration and the DDW state, the phase transition between 
them can been seen most directly from the particle num- 
ber distribution in the real place, as shown in Fig|6l 
where we set t = Q.\,K = 1.0, L^, = 26. In the phase 
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FIG. 6: The particle distribution in the real space in a 2 x 
26 ladder, we set t=0.1, K=1.0 and (a)V=0.01 (b)V=0.41, 
(c)V=0.44, (d)V=0.46, (e)V=3.0 



separation region (V=0.05), due to the open boundary 
condition, the boson prefer to get together in the center 
of the ladder to form a half- filling DBL and the rest part 
is empty^. Deep in the DDW state (V=3.0), there is a 
three-period crystal structure of the particle number dis- 
tribution. There is a quantum phase transition between 
these two phases. 

To study the properties of the quantum phase transi- 
tion, we introduce AA^ — p{Lx/2) — p{Lx/2 — 1) as the 
order parameter to characterize the density wave state, 
where p{i) has been defined above. The dependence of 
AiV on the NN interaction V is shown in FiglTJa), where 
we set t = 0.1 and K = 1.0 for simplicity. In a perfect 
DDW state {V — oo), AA^ = 0.5, while in the opposite 
limit (V=0), a phase separation means AiV = 0, (in our 
case K/t=10, thus the ground state in K=0 is phase sep- 
aration rather than DBL^^). FiglTl^a) indicate a strong 
signature of a continuous quantum phase transition, in- 
stead of a simple first-order phase transition caused by 
energy level crossing. To verify this point, we also calcu- 
late the dependence of the average ground state energy 
per boson (Fig. El^b)), its first as well as the second or- 
der derivative (FiglTKc) and (d)) on the NN interaction V. 
We notice that there is no discontinuity in the first order 



derivative (§^) while a sharp peak appears in the second 
order derivative of the ground state energy (^^), which 



indicates that a second order phase transition occurs at 
the point K = 0.43 ~ 0.44. 

The pairing between two bosons have recently at- 
tracted considerable attentions, while most of previ- 
ous pairing mechanisms are based on the direct at- 
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FIG. 7: (a)The dependence of AA'^ on V in a two-leg ladder, 
we set t=0.1, K=1.0; (b)The average energy per boson E .vs. 

V; (c) .vs. V; (d)|f .vs. V. 



tractive interspecies interactions tuned by Feshbach 
Resonance^^^, or on a three-body onsite hardcore 
constraint'^''. All of these mechanisms are due to the 
uniqueness of the cold atomic system and have no coun- 
terpart in traditional condensed matter physics. In this 
paper, we proposed a novel pairing mechanism for bosons 
via strong four-site ring exchange interaction, which is 
also thanks to the unique feature of the ultracold atoms in 
optical lattice. Recently a proposal has been provided to 
experimentally detect these bosonic pair—, which would 
be helpful to detect the bosonic pair in our case. 

In summary, we propose a strongly correlated bosonic 
Hamiltonian with four-site ring exchange interaction. We 
focus on the parameter region V ^ K ^ \t\ > and in- 
vestigate the exotic phases with d-wave symmetry emerg- 
ing at different filling factors. A physical realization of 
the Hamiltonian as well as the parameter region has also 
been discussed. 

The authors are grateful to Congjun Wu for helpful dis- 
cussion. ZC is supported by AROW911NF0810291. The 
work is also supported in part by NSF-China, MOST- 
China and US-DOE-FG02-04ER46124 (XCX). 
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